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ABSTRACT 

Raman scattering by atomic hydrogen converts the UV continuum around Ly/3 into 
optical continuum around Hct, and the basic atomic physics has been discussed in 
several works on symbiotic stars. We propose that the same process may operate in 
active galactic nuclei (AGN) and calculate the linear polarization of the broad emission 
lines Raman-scattered by a high column neutral hydrogen compnent. The conversion 
efficiency of the Raman scattering process is discussed and the expected scattered 
flux is computed using the spectral energy distribution of an AGN given by a typical 
power law. The high column Hi component in AGN is suggested by many observations 
encompassing radio through UV and X-ray ranges. 

When the neutral hydrogen component with a column density ~ 10^^ cm“^ is 
present around the active nucleus, it is found that the scattered Ha is characterized 
by a very broad width ~ 20, 000 km s“^ and that the strength of the polarized flux is 
comparable to that of the electron-scattered flux expected from a conventional unified 
model of narrow line AGN. The width of the scattered flux is mainly determined 
by the column density of the neutral scatterers where the total scattering optical 
depth becomes of order unity. The asymmetry in the Raman scattering cross section 
around Ly/3 introduces red asymmetric polarized profiles around Ha. The effects of 
the blended Ly/3 and O vi 1034 doublet are also investigated. 

We briefly discuss the spectropolarimetric observations performed on the Seyfert 
galaxy IRAS 110548-1131 and the narrow line radio galaxy Cyg A. Several predictions 
regarding the scattering by a high column neutral hydrogen component in AGN are 
discussed. 
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1 INTRODUCTION 

The spectra of active galactic nuclei (hereafter AGN) are 
characterized by many prominent emission lines. The width 
of the emission lines is one of the important parameters to 
classify these active galaxies. For example, type 1 Seyfert 
galaxies possess both broad and narrow emission lines and 
type 2 Seyfert galaxies exhibit only narrow emission lines 
(e.g. Osterbrock 1989). Radio galaxies are also classified into 
broad line radio galaxies and narrow line radio galaxies in a 
similar way. The width of the broad emission lines (BEL’s) 
often ranges up to several thousand km s“^ and the narrow 
emission lines (NEL’s) show an order of magnitude smaller 
width in both type Seyfert galaxies. 

A unification model of Seyfert galaxies suggests that 
Seyfert 2 nuclei also have the broad emission line region 
(BELR), which is blocked from the observer’s line of sight. 
In this unification scheme, the BEL photons can be scattered 
from other lines of sight to enter the observer’s line of sight. 
Because the scattered component is expected to be highly 


polarized, spectropolarimetry plays an important role to test 
the unification models (e.g. Antonucci 1993). 

A similar behavior is seen in broad absorption line 
quasars, in which the broad absorption trough is partially 
filled by photons scattered from other lines of sight. These 
photons are believed to be scattered resonantly by a fast 
moving medium just outside the BELR, and several spec¬ 
tropolarimetric observations using large telescopes such as 
the Keck 10 m telescope have been performed (e.g. Cohen 
et al. 1995, Goodrich & Miller 1995, Lee & Blandford 1997). 

Spectropolarimetry of Seyfert 2 galaxies performed by 
several researchers confirmed that the BEL’s are shown in 
the polarized flux as expected from the unification scheme 
(e.g. Antonucci & Miller 1985, Miller & Goodrich 1990, An¬ 
tonucci 1993). The constant position angle and degree of 
polarization favor strongly the model, in which the electron 
plasma located in high latitude scatters the radiation origi¬ 
nated from the BELR into the observer’s line of sight. 

Tran (1995 a,b,c) presented the spectropolarimetric 
data on a number of Seyfert 2 galaxies, which exhibit broad 
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emission lines in the polarized flux. He noted that the 
Balmer lines are more highly polarized than the adjacent 
continuum even after all the corrections due to the star light 
and the interstellar medium. He proposed an explanation by 
assuming that there is still unpolarized continuum having a 
flat spectra not included in the correction procedures. This 
component provides more dilution to the continuum leaving 
the lines more polarized. However, higher stellar contribu¬ 
tion as suggested by Koski (1978) would result in wavelength 
independent polarization (Antonucci, 1993). 

Another interesting and important spectropolarimet- 
ric observation is presented by Ogle et al. (1997), who ob¬ 
served the narrow line radio galaxy Cyg A. They found that 
the Ha polarized flux has a very broad width exceeding 
20000 km s“^. Young et al. (1993) also found very broad 
polarized Ha from their spectropolarimetry of the Seyfert 
2 galaxy IRAS 110548-1131, where the full width at zero 
intensity is 16800 km s“^. 

The X-ray observations of Seyfert galaxies have been 
performed to elucidate the physical properties of the highly 
ionized component plausibly related with the scattering 
agents of the Seyfert 2 galaxies. The warm absorbers having 
a column density of Nh ~ 10^^ — 10^^ cm~^ are invoked 
to explain the soft X-ray observations of Seyfert galaxies. 
On the other hand, the existence of the neutral component 
with a similar column density is controversial (e.g. Marthur 
1997). 

Conway & Blanco (1995) used the VLA to obtain the H l 
column density N[fi ~ iq22-23 toward the compact 

nucleus of Cyg A, by measuring the 21 cm absorption. The 
neutral hydrogen is concentrated near the nucleus within 
several parsecs and they suggested that atomic hydrogen is 
an important constituent of the molecular torus obscuring 
the BELR of the narrow line radio galaxy. Their result is 
compatible with an X-ray spectroscopy by the Ginga satel¬ 
lite (e.g. Ueno et al. 1994) suggesting a highly optically thick 
component around the compact nucleus. 

Strongly polarized broad emission features around 
AA6830, 7088 in a fair fraction of symbiotic stars were re¬ 
cently identified to be Raman-scattered O vi AA1032, 1038 
by atomic hydrogen of large column density (e.g. Schmid 
1989, Nussbaumer et al. 1989). The absorbed UV photons 
with wavelength Ai < Aa = 1216 A by a hydrogen atom 
in the ground state Is can be re-emitted either coherently 
with wavelength A/ = Ai (the Rayleigh scattering) or in¬ 
coherently with A/ = (A~^ — Aa^)“^ (the Raman scatter¬ 
ing). Therefore, UV continuum photons near Ly/1 are scat¬ 
tered and transformed into optical photons around Ha by 
the Raman process. One of the most important characteris¬ 
tics of the Raman scattering process is the enhancement of 
the Doppler shift so that the scattered waves appear to be 
shifted by a factor of ~ A//Ai, and this leads to a forma¬ 
tion of a very broad line feature from a moderately broad 
incident radiation source. 

The Raman process can be important in the systems 
having both a strong UV incident radiation source and a 
neutral hydrogen component of large column density. The 
optimum conditions for the Raman scattering may be pro¬ 
vided by AGNs, where the spectral energy distribution 
(SED) is characterized by the big blue bump in the UV range 
and neutral components with a high column density may be 
found near the broad line region. As in symbiotic stars, the 


Raman scattering may produce observable features around 
Ha both in the profile and in the polarization. 

In this paper we investigate the possible observational 
features expected if the Raman process operate in AGN. 
The paper is composed as follows. In section 2, we give a 
model description and make an estimate of the strength of 
the Raman-scattered flux using a standard AGN model. In 
section 3, the Monte Garlo result is presented and compared 
with the electron-scattering model. We discuss the applica¬ 
tion of the Monte Garlo result to a couple of narrow line 
AGN and the observational implications in section 4. We 
summarize in the final section. 


2 MODEL 

2.1 Scattering by a high column neutral hydrogen 

The polarized fluxes of many Seyfert 2 galaxies are charac¬ 
terized by a constant position angle and a constant degree of 
polarization, which strongly indicate that they have the elec¬ 
tron scattering origin. The width of the broad polarized flux 
puts a constraint on the electron temperature Te 10® K. 
Assuming a covering factor ^0.1 and a similar luminosity of 
the hidden broad line region of Seyfert 2 galaxies to that of 
Seyfert 1 galaxies, one can deduce that the electron column 
density is of order 10^^ cm“^ (e.g. Miller & Goodrich 90). 

Recent X-ray observations on several Seyfert galaxies 
have shown that there exists a highly ionized component of 
column density 10^^ —10^^ cm“^. The warm absorber model 
has been introduced to account for the X-ray absorption as 
well as the UV absorption properties of the Seyfert galaxies 
(e.g. Marthur 1997, Kriss et al. 1995). Therefore the warm 
absorber may provide a good candidate for the free electron 
scatterers of the UV photons originating from the broad line 
region. 

A high column neutral component has been suggested 
from X-ray observations of Seyfert galaxies (Turner et al. 
1996, Nandra et al. 1994, Netzer 1993). The physical nature 
of the X-ray absorbing material in AGN is still controversial, 
and the possibility that both cold and warm absorbers may 
co-exist is not excluded. 

Many X-ray observations of Seyfert 2 galaxies show that 
Gompton reflected components seen in the spectra in the 
band 0.1 keV - 10 keV are compatible with the assumption 
of the molecular torus supplying the cold optically thick ma¬ 
terial. It is not certain whether the neutral absorbers sug¬ 
gested for the X-ray observations of Seyfert 1 galaxies are 
also originated from the molecular torus or constitute an 
independent component in the center part of AGN. 

More direct information about the amount of neutral 
hydrogen can be obtained from radio observations using the 
21 cm line from the hyperfine structure of hydrogen atom. 
Dickey (1986) investigated the content of neutral hydrogen 
in 19 active spiral galaxies and obtained values in the range 
Nfji ~ 10^°“^^ (Ts/lOO K) cm“^. He also mentioned the 
possibility that the spin temperature Ts may reside in a 
large range 10^“^ K, in which case a column density as high 
as 10^® cm“^ is also allowed. 

Similar observations have been made by a number of 
researchers including Gallimore et al. (1994), Mundell et 
al. (1995). It seems that the results of these observations 
produce more or less similar values for the neutral column 



Polarization of Ha in Active Galactic Nuclei 3 


density around the galactic nuclei. 

In this subsection, we investigate the scattering prop¬ 
erties of atomic hydrogen. We will assume a neutral com¬ 
ponent along the normal direction to the accretion disk 
with a covering factor of about 0.1 and column densities 
10^1-23 cm-2. 

Interactions of electromagnetic waves with atoms are 
described by the Kramers-Heisenberg formula (e.g. Sakurai 
1967). The scattering types are divided into the Rayleigh 
scattering and the Raman scattering, where in the former 
case the initial and the final states of the scatterer are the 
same and the scattering is coherent. On the other hand, in 
the case of the Raman scattering the final state of the scat¬ 
terer differs from the initial state and therefore the scat¬ 
tered photon has different frequency from that of the inci¬ 
dent one. The explicit numerical computation of the relevant 
scattering cross sections has been done by many researchers 
(Gavrila 1967, Saslow and Mills 1969, Schmid 1989, Isliker 
et al. 1989, Sadeghpour & Dalgarno 1992, Lee & Lee 1997a). 

Following Lee & Lee (1997a), the Rayleigh and the Ra¬ 
man scattering cross sections are given by 



where the subscript v = {Ray, Ram} represents the scat- 
2 

tering type, ro = ^ ^2 is the classical electron radius, uj 
and e are the angular frequency and the polarization vector 
for the incident wave and the primed quantities correspond 
to the scattered wave. The matrix elements for the bound 
states are 
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where the double-bar (or reduced) matrix elements for the 
position operator are found in Berestetskii et al. (1972) (see 
also Bethe 1967). Here, u!ni{i = 1, 2} is the transition angu¬ 
lar frequency between the level n and Is or 2s depending on 
i. The contributions to the matrix elements Mfay, Hfam 
from the continuum states are similarly given (e.g. Saslow 
& Mills 1969). 

In Fig. 1, we show the total (Rayleigh-fRaman) and 
the Raman scattering optical depths for a neutral hydrogen 
slab of column density Nr \ = 10^^ cm“^ as a function of the 
incident wavelength. We present the incident wavelength Ai 
by the velocity deviation defined by AV = c(Ai — XfPj/^p- A 
close inspection reveals that the total scattering cross section 
is also slightly asymmetric. More quantitatively, in the range 
—850 km s“^ < AH < 1200 km s”*^, the total scattering 
optical depth exceeds unity. 

The ratio of the Raman scattering cross section to that 
of the total scattering cross section is also plotted in Fig. 1. 


It is particularly notable that the ratio is asymmetrical. The 
ratio increases from 0.03 to 0.3 in the wavelength interval 
1010 A < A < 1040 A. This implies that the conversion 
efficiency from UV around Ly/3 to optical around Ha is much 
larger in the red part than in the blue part. 

From Fig. 1, we see that for a slab of column Nri = 
10^® cm“^ incident photons with the velocity difference sat¬ 
isfying —3300 km s“^ < AH < 3100 km s“^ (or 1014 A < 
A < 1036 A) have total scattering optical depth Ttot > 1. The 
wavelength range considered in this work is described by the 
Lorentzian wing part, and the Doppler shift corresponding 
to the total scattering optical depth of unity depends ap¬ 
proximately on the square root of the column density, i.e., 

|AH| ~ 2000 km s"\ (2.4) 

where N 22 = Nri/{1{P^ cm“^). 

The incoherent nature of the Raman scattering process 
introduces an enhancement of the Doppler shift by the rela¬ 
tion 



(e.g. Schmid 1989). In the case of incident photons hav¬ 
ing wavelengths similar to that of Ly/3, the wavelength ra¬ 
tio A//Ai ~ 6.4, and therefore the incident photons hav¬ 
ing a width ~ 2000 km s“^ around Ly/3 will result in a 
broad feature around Ha with a width ~ 13,000 km s“^. 
This implies that a broad line-like scattered feature with 
AH ^ 15, 000 km s“^ can be formed from a neutral hydro¬ 
gen component with Nri ^ 10^^ cm~^. 

In this work, we do not consider the contribution of 
the Ha line photons scattered by hydrogen atoms in excited 
states n — 2, because the velocity width for this process 
to be important is a few times the thermal velocity width 
depending on the population of the excited state hydrogen 
atoms and is much smaller than the typical velocity scale 
of the BELR. Therefore, this process has a very small band 
width compared with that of the Raman scattering process. 
However, if we have a neutral scattering component cover¬ 
ing a large velocity space ranging ^ 1000 km s“^, then this 
process should not be neglected. 

2.2 Spectral Energy Distribution 

In order to assess the flux Raman-scattered by neutral hy¬ 
drogen and compare it with the flux scattered by free elec¬ 
trons, it is essential to know the SED of the AGN both in 
UV and in optical ranges. A typical AGN spectrum is char¬ 
acterized by the so called big blue bump in the ultraviolet 
range and the specific luminosity shows a peak at ~ 10^® Hz 
(e.g. Mathews & Ferland 1987). It is conventional to repre¬ 
sent the SED as a power law with the power index a so that 
we have 



Zheng et al. (1997) investigated more than 200 quasar 
spectra obtained from the Hubble Space Telescope (HST) 
and constructed composite quasar spectra. They concluded 
that the power law index a = —0.99±0.05 in the continuum 
between 1050 A and 2200 A. 

Natali et al. (1998) also investigated the optical-UV 
continuum SED using the spectra of 62 QSOs emphasizing 
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Figure 1. 

Total (Rayleigh+Raman) and Raman scattering optical depths for a neutral hydrogen slab of column density Nhi = 10^^ cm“^. 
The horizontal axis shows the incident photon wavelength in units of A. The solid line shows the total scattering optical depth and the 
Raman scattering optical depth is represented by the dotted line. The dashed line represent the ratio of the Raman scattering cross 
section to that of the total scattering cross section. The points marked by a cross are the velocity shifts having the total scattering optical 
depth of unity for Nhi = 10^^ cm“^, and the points marked by a circle are for Nhi = 10^^ cm“^. 


that the local nature of the power-law approximation. They 
obtained the mean values of the power-law index change 
from a ~ 0.15 at A > 3000 A to a ~ —0.65 at A < 3000 A. 

We assume that the incident spectrum in the UV and 
the optical ranges is given by a double power law, i.e., 


if ^>*^0 

if 


(2.7) 


where uo = 10^® Hz, and we set a\ = —1, and 02 = 0.15. 

According to the prescription Eq. (2.7), the incident 
photon number flux Oa is related with the speciflc luminos¬ 
ity 
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and therefore, the continuum number flux near Ha is given 
by 


4>a(A ~ \Ha) ~ 0.4 X 4>a(A ~ ^Lp) 


(2.9) 


Lee & Lee (1997b) discussed the basic properties of the 
Raman scattering process, and gave an empirical formula of 
the efficiency eRam for the conversion of the UV photons into 
the optical photons. They showed that the efficiency is quite 
large when the total scattering optical depth is moderately 
large, or, Tscat ^ 1. Furthermore, the conversion efficiency is 
sensitively dependent on the ratio of the Raman scattering 
cross section to the total scattering cross section. 

Using the relation 



( 2 . 10 ) 


and combining it with Eq. (2.9), we obtain the Raman- 
scattered continuum flux at around Ha 

4>f“™dA/ ~ CNtRam^X,d\i ^ 

~ 0.06 CNtRam^XfdXf. 

Here, Ujv is the covering factor of the neutral component 
having the total scattering optical depth ^ 1, and A/ ~ 
6563 A , Ai ~ 1025 A , respectively. 

The continuum number flux $Th that is Thomson scat¬ 
tered near Ha depends on the electron column density and 
the covering factor. Hence, assuming that the scatterers have 
a small Thomson scattering optical depth tw < 1, we may 
write 


4>r ~CwTm4>A^, (2.12) 

where Cw is the covering factor of the ionized scattering 
medium. Here, we ignore the line broadening by the thermal 
motions of the electrons. 

Therefore, the Raman-scattered continuum flux around 
Ha is compared with the Thomson-scattered continuum flux 
by the relation 

Under the assumption that the electron plasma has a 
column density 10^^ cm“^ (or tw = 0.6 x 10“^), and that 
the neutral hydrogen has the column of Nri = 10^^ cm~^, 
the continuum flux that is Raman-scattered from Ly/3 with 
a width of 2000 km s“^ is comparable to the Thomson- 
scattered at around Ha with a width of 15,000 km s“^. 
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Because enam is a sensitive function of the column density 
and the frequency, the Raman-scattering material can be 
regarded as a scattering-mirror with a narrow band width 
centered at the resonance frequencies (e.g. Lee & Lee 1997 
a,b). 

In an AGN spectrum, both Ly/3 and Ha are very strong 
lines and especially Ly/3 is often blended with O vi 1034 
to form a very broad feature around 1030 A and thus the 
line photons are also an important constituent of the UV 
incident radiation (e.g. Laor et al. 1994, Laor et al. 1995, 
Netzer 1990). 


2.3 Model and the Monte Carlo Method 


We briefly describe the model and the Monte Carlo code 
to compute the polarization due to the neutral component. 
In this paper the neutral scatterers are placed along the 
symmetry axis for simplicity and ease of a direct comparison 
with electron-scattered flux. A more realistic model may be 
suggested by placing the neutral scatterers in the molecular 
torus, the top part of which is accessible by an observer 
with low latitude. If the molecular torus is approximated by 
a perfect cylinder with an inner radius a and a height h (e.g. 
Pier & Krolik 1992), then the total solid angle subtended by 
the accessible part of the torus by an observer with the line 
of sight ko with respect to the light source at the origin is 
given by 


If 


dz d<f) 


z>|+2ocoteo cos,/> [a2 -I-2;2]3/2 

|.7r/2 

= 4 cot do du u sin u 


(h 

1 -I- ( ^-2cot 9o cos u\ 


(2.14) 


Here, we take the cylinder axis as 2 ;—axis and 60 = cos“^ ko • 
z is the angle between the line of sight and a—axis. In 
Fig. 2a, a schematic geometry is shown and the hatched 
region marked by ‘A’ is the accessible part by the observer. 

For parameters h/2a = 1 and 6 q — cos“^ 0.5, a sim¬ 
ple numerical integration of Eq. (2.14) gives the solid angle 
AH/47r = 0.08. Hence, this particular geometry mimics a 
neutral cloud located at high latitude to the observer with 
the covering fraction ~ 10 percent. However, a tapered disk 
model investigated by Efstathiou et al. (1995) would give 
much smaller solid angle. Therefore, only after a detailed 
geometry of the cicumnuclear region is specihed with the 
appropriate photoionization computation, the covering frac¬ 
tion of the neutral scatterers may be assessed. In this paper 
we choose a very simple geometry as depicted in Fig. 2b to 
hnd the main features of Raman scattering by neutral hy¬ 
drogen without invoking many free parameters necessary to 
specify a detailed geometry. In this geometry, the position 
angle of the polarized radition is perpendicular to the sym¬ 
metry axis, which is regarded as the radio axis of the active 
nucleus. 

As shown in Fig. 2b, the light source is located at the 
center and blocked by an opaque component lying equatori- 
ally, and the neutral hydrogen component lies alogn 2 —axis 
that is coincident with the normal direction of the accre¬ 
tion disk. The geometry of the scatterer is assumed to be an 


oblate spheroid with the minor axis aligned to the 2 —axis. 
Along the major axes the column density is assumed to be 
2 Nh I and inside the scattering region, the number density is 
assumed to be hxed, which enables us to measure the physi¬ 
cal distance in terms of column density. The covering factor 
Cn of the neutral scatterers is set to be 0.15. 

We prepare the UV incident continuum in accordance 
with the double power law given by Eq. (2.7). We also pre¬ 
pare the line photons around 1030 A corresponding to the 
blended Ly/3 and O vi. We combine the three gaussians cen¬ 
tered at 1025 A, 1032 A and 1038 A with the same width 
of 3000 km s“^ and strength ratio of 1 : 1.5 : 1.5. The total 
blended line flux is renormalized so that it has an equiva¬ 
lent width of 30 A (e.g. Korista et al. 1997). The possibility 
for intrinsic anisotropy in the radiation field in the BELR 
has been suggested and observational evidence is presented 
for NGC 1068 (Miller et al. 1991, Hough and Young 1995). 
However for simplicity it is assumed that the light source is 
isotropic. 

As a UV photon reaches the scattering region, it may 
be first Rayleigh-scattered several times before it is Raman- 
scattered. It is assumed that the continuum opacity near 
Ha is negligible and therefore any Raman-scattered photon 
leaves the scattering region. Here, the thermal motion of the 
scatterers is not important because the bulk velocity of the 
BELR is much larger. 

The polarization of the line photon resonantly scattered 
by a hydrogen atom is dependent on the incident wavelength 
in the rest frame of the scatterer. However, most UV pho¬ 
tons are scattered in the damping wing regime where the 
deviation from the resonance transition frequency is suffi¬ 
ciently large to ensure that the scattering is characterized 
by the classical Rayleigh phase function. This is a direct con¬ 
sequence that the fine structure splittings in the np states 
are negligible compared with the frequency deviation cor¬ 
responding to the damping wing regime (e.g. Stenflo 1980). 
Because the final states for both the Rayleigh scattering and 
the Raman scattering are s states, the same Rayleigh phase 
function is used to obtain the wave vector and the polar¬ 
ization associated with the scattered photon for both the 
Rayleigh and the Raman scatterings. 


3 RESULTS 

3.1 Continuum and Basic Properties 

In Eig. 3, we show the scattered flux, the polarized flux 
and the linear degree of polarization around Ha resulting 
from Raman scattering by neutral hydrogen of column den¬ 
sities in the range Nfji — ® cm~^. The UV incident 

spectrum consists of only the continuum given by Eq. (2.7), 
and we collect the photons with wave vector k satisfying 
0.45 < k • z < 0.55, where 2 —axis is the symmetry axis of 
the AGN. 

The polarized flux appears to be a normal Ha line with 

FWHM ~ 1.6 X 10"^ cm~^)^/^ km s"\ (3.1) 

This corresponds to a few times the width shown in 
Eq. (2.4). The scattered flux shows a peak at the line 
center, and has a red excess that becomes more conspic¬ 
uous as the column density increases. The polarized flux 
is broader than the scattered flux and shows a red-shifted 
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Figure 2a. 

A schematic diagram for a cylindrical molecular torus (see Pier Sz Krolik 1992). The inner radius is a and the height is h. 2 —axis is 
chosen to coincide with the symmetry axis and the central continuum source is located at the origin. The hatched region marked by ‘A’ 
can be accessible by the observer with the line of sight ko- Assuming the the observer is in the x — z plane, the region ‘A’ is described 
by the relation z > 2a cot Oq cos where Oq = cos“^ ko ■ z. 


Figure 2b. 

A schematic diagram for the scattering model. The continuum light source lies in the center and the broad line photons are generated 
in the vicinity of the continuum source. The neutral hydrogen component is located along the 2 :—axis, which is coincident with the normal 
direction to the accretion disk. 
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Figure 3. 

The scattered flux, the polarized flux and the linear degree of polarization around Ho Raman-scattered by a neutral hydrogen scattering 
component of column densities N[fi = 10^^, 10^^ ®, 10^^, 10^^ ® cm“^. The UV incident spectrum is the continuum given by Eq. (2.7), 
the photons with wave vector k satisfying 0.45 < k • z < 0.55 are collected, where z axis is the symmetry axis of the AGN. The solid line 
represents the scattered flux the dotted line stands for the polarized flux multiplied by 6 for comparison with the scattered flux and the 
dotted line with error bars for the degree of polarization. The error bars are l-cr obtained from the Monte Carlo computation. 


peaks. The locations of the peak of the polarized flux are 
given by AVpeak ~ 5600 km s“^ for Njji = 10^^ ® cm“^, 
and AVpeak ~ 1700 km s“^, and ~ 800 km s“^ for N^i = 
10^^ cm“^ and 10^^'® km s“^, respectively. 

Near the line center, the scattering optical depth is very 
large, and therefore a large scattering number is needed to 
escape from the scattering region. As the scattering number 
gets larger, the conversion from Ly/3 to Ha becomes more 
efficient, and hence we obtain a stronger Raman flux. 

Furthermore, the conversion efficiency is dependent on 
the ratio of the Raman scattering cross section to that of 
the total scattering cross section in a highly sensitive way. 
Therefore, as the column density becomes large, the mean 
scattering number also increases and the strength of the 
Raman scattered flux is enhanced preferentially to the red 
part where the ratio of the Raman scattering cross section 
is larger as shown in Fig. 1. 

An illustrative explanation is given by Lee & Lee (1997 
b) in their analysis of the Raman scattered flux in symbiotic 
stars. It is notable that the red asymmetric polarized profile 
in the high column case is originated from the atomic physics 
rather than extrinsic factors such as kinematics. 

On the other hand, the largest linear degree of polariza¬ 
tion is usually obtained when the mean scattering number 
is ~ 1, because a large number of scatterings randomize the 
radiation field associated with the scattered photons and 
no scattering yields any polarization. However, in the case 
of the Raman scattering, where the scattered photons are 
not mixed with the unscattered flux due to the incoherence, 
the degree of polarization monotonically decreases as the 


mean scattering number increases. Therefore, the polariza¬ 
tion tends to decrease as the total scattering cross section 
increases. 

Fig. 3 shows that irrespective of the column density, 
the degree of polarization is the smallest at the Ha cen¬ 
ter, where the total scattering cross section is the largest. 
It is also interesting that the minimum degree of polariza¬ 
tion at the center is approximately constant independent of 
the scattering column density. This behavior is explained 
by noting that the Raman flux is mainly contributed from 
the first few scatterings and that negligibly polarized Ra¬ 
man photons contribute little because of low probability of 
surviving a large number of Rayleigh scatterings. Near the 
center, the total scattering optical depth is large enough 
even for Npji — 10^^ cm“^ that the polarized flux and the 
degree of polarization become saturated. 

Because the total scattering optical depth rapidly de¬ 
creases as the wavelength deviation from the center wave¬ 
length gets larger, the scattered flux becomes negligible. 
This is also illustrated by the large error bars in the figure. 
It is interesting to note that the polarized flux, the total 
scattered flux and the linear degree of polarization are de¬ 
pendent on the wavelength and the column density of the 
scatterer, as may be manifested by fine structures in the 
spectropolarimetric observations. 

3.2 Line Contribution and Electron Scattering 

In Fig. 4 are shown the same quantities as in Fig. 3 
except for the UV incident radiation that includes the line 
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Figure 4. 

The same quantities as in Fig. 3 except for the UV incident radiation includes the line contribution from the Ly/9 and O VI doublet. 


contribution from the Ly/3 and O vi doublet. A gaussian 
profile for each line is assumed. With a continuum flux Fc 
and a line flux 

Fi = Fo exp[-(A - Ao)V2<rf] (3.2) 

the equivalent width EW is given by EW = -y/TTailEo/Fc)- 
We take ai = Ao(3000 km s“^/c) for each line. The equiv¬ 
alent widths of the components of the blended feature are 
given by EWhyp = 7.5 A, EWqyi = 22.5 A, where the dou¬ 
blet O VI is assumed to be contributed equally from A1032 
and A1038 components. 

The qualitative features shown in Fig. 4 are similar to 
those in Fig. 3. The main difference is the strength of the po¬ 
larized flux, due to the increase of the incident UV photons 
by a factor of ~ 2. Because the Raman-scattering mirror 
has a narrow band width, the overall variation of the inci¬ 
dent UV flux over this mirror-width is small and therefore 
we obtain qualitatively similar polarized fluxes. 

However, because the blended line feature is centered 
near 1030 A, the red side of Ly/l is stronger than the blue 
side by a small amount. This leads to more red-asymmetric 
Raman-scattered flux, which is barely apparent in the fig¬ 
ure. The small effect of O vi is due to the small scattering 
cross section, and the major contribution to the polarized 
flux is provided by the Ly/3 line photons which possess a 
large scattering cross section. When N^i ^ 10^^ cm“^, the 
Raman scattering is not important and the asymmetry in 
the profile is negligible. 

In Fig. 5 we compare the Raman-scattered polarized 
flux with the polarized flux scattered by an electron plasma 
located in the same position with the same covering fac¬ 
tor as the neutral hydrogen scatterer. The column density 
Ne of the free electron is assumed to be Ne = 10^^ cm“^. 
The incident flux consists of the continuum described in 


Eq. (2.7) and the broad Ha given by a gaussian with 
(Ji = 3000 km s“^ and the equivalent width EWhc ~ 100 A. 

If we assume that the thermal velocity of the electron 
plasma is much smaller than the width of the broad Ly/?, 
then the line broadening by thermal electrons is neglected. If 
the Thomson scattering optical depth is also small, then the 
Thomson scattered flux is mainly contributed from singly 
scattered photons. Under these assumptions, the electron- 
scattered flux is characterized by a constant degree of po¬ 
larization and therefore, the polarized flux has almost the 
same shape as the that of the scattered flux 

Therefore, in this optically thin limit, the degree of po¬ 
larization pTh of the Thomson scattered flux is given by 


PTh 


1 

27r(l — He) 


dp, d(f> 


Pi - P2 
Pi +P2 


(3.3) 


where pi = (cos 9o sin 0)^-1- (cos 6 o cos 8 cos <j) + sin 9o sin 9)'^ 
and p 2 = cos^ (p + cos^ 9 sin^ p are the scattered compo¬ 
nents that are polarized parallelly and perpendicularly to 
the plane spanned by the observer’s line of sight and the 
symmetry axis (a—axis) of the system, respectively (e.g. 
Brown & McLean 1977). Here, 9o is the angle between the 
2 axis and the observer’s line of sight, p^ is the cosine of the 
half-angle of the conical region which the electron plasma 
subtends with respect to the source. 

For 9o ~ cos“^ 0.5 and p^, = 0.7, a simple numerical 
integration gives pTh = 0.39, which agrees well with the 
largest degree of polarization obtained for the Raman scat¬ 
tered flux shown in Fig. 3. 

The electron-scattered flux around Ha is comparable 
with the Raman scattered flux for Nui = 10^^ ® cm“^. 
We include an unpolarized continuum of strength <l>„c = 
5 X 10“®<1>A as suggested by Tran (1995 c) to compute the 
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synthetic linear degree of polarization represented by the 
long-dashed line in Fig. 5. With the polarization-diluting 
continunm, the degree of polarization shows a bump-like fea¬ 
ture around Ha, observed in a number of Seyfert 2 galaxies 
(e.g. Tran 1995 a,b,c). It is particularly noted that the polar¬ 
ization bumps also have a dip in the center which is ascribed 
to the decrease of polarization arising from a large number 
of Rayleigh-Raman scatterings. 


4 OBSERVATIONAL IMPLICATIONS 

As mentioned in section 2.2, the Raman-scattering material 
may be regarded as a scattering-mirror with a very narrow 
band width centered at the resonance frequencies and fur¬ 
thermore the band width provides a natural measure of the 
column density of the neutral hydrogen. Therefore, when the 
incident flat continuum is reflected by the Raman-scattering 
mirror, the scattered continuum appears a line-like feature 
around the transition frequency Vi — VLya- 

Stark et al. (1992) presented a Hi survey of our galaxy 
using a beam with a size of FWHM= 2° and showed that 
the HI column density toward the galactic center amounts to 
10^^ cm“^. Gallimore et al. (1997) presented the Hi column 
density toward the radio jet of Mrk 6, and showed that neu¬ 
tral hydrogen of column ~ 10^^ cm“^ is present around the 
active nucleus. X-ray observations of many AGN reveal that 
X-ray absorbing material of A ~ ]^q 20-23 around the 

nuclear region and the physical nature, the exact location 
and kinematics are currently being investigated intensively. 

Ogle et al. (1997) performed a spectropolarimetry of 
the narrow line radio galaxy Cyg A and found that the 
full width at half maximum of the polarized Ha flux is 
26000 km s“^. They noted that the broad emission lines of 
some radio loud active galaxies show a large velocity width 
up to 30000 km s“^ with double-peaked profiles (e.g. Era- 
cleous & Halpern 1994). However, Antonucci et al. (1994) 
observed the same object and measured the width of the 
broad emission line Mg H ~ 7500 km s“^, which is 3.5 times 
smaller than that of the polarized Ha. 

The thermal width of an electron plasma with the tem¬ 
perature Te is given by 

~ 400 km s-\ (4.1) 

where Te 4 = Te/(10^ K). If the width of the Ha in the polar¬ 
ized flux is interpreted to be caused by hot electrons, then 
the temperature of the electron plasma Te ~ 5 x lO’^ K (e.g. 
Weymann 1970). And there may also exist other polarized 
lines with a similar width, and it is difficult to explain the 
width of Mg H line observed by Antonucci et al. (1994). 

On the other hand, if the Raman-scattering processes 
by the neutral component operate, then it may provide an 
interesting interpretation, i.e., the large width may not be 
associated with the Doppler effect but mainly an effect of 
the atomic physics associated with the Raman scattering. 

Conway & Blanco (1995) measured the column den¬ 
sity of HI toward the compact radio nucleus of Cyg A and 
obtained a value > 2.54 ± 0.44 x 10^® Ts cm~^, where the 
spin temperature Ts may go up as high as 1000 — 8000 K de¬ 
pending on the physical state of the putative molecular torus 
obscuring the compact nucleus. Therefore the existence of a 
high column neutral hydrogen with Nfn ^ 10^^ cm“^ is not 


ruled out. The Raman-scattering by this component may in¬ 
troduce a broad polarized flux with width ~ 2 x 10^ km s“^ 
redshifted by ~ 2000 km s“^ from incident UV radiation 
around Ly/I with width > 3000 km s“^. Therefore, if we 
take the width of Mg H 2800 as a representative velocity 
scale of the BELR, then the spectropolarimetry of Cyg A 
by Ogle et al. (1997) is consistent with this interpretation. 

Young et al. (1993) presented their spectropolarimetric 
observation of the Seyfert 2 galaxy IRAS 110548-1131 (see 
also Young et al. 1994). They found a very broad Ha line 
in the polarized flux with a full width at zero intensity of 
16800 km s“^ and FWHM= 7600 km s“^, which, they note, 
is wider than that of most Seyfert 1 galaxies. Another in¬ 
teresting point is that the peak of the broad polarized Ha 
is red-shifted by 900 ± 400 km s“^ with respect to the nar¬ 
row Ha emission, which they ascribe to an outflow of the 
scatterers. 

If we apply our result, the most plausible column den¬ 
sity for IRAS 110548-1131 is ~ 10^^'® cm“^, which gives the 
FWHM of the polarized flux ~ lO'^ km s“^ and the redshift 
of the peak ~ 800 km s“^. 

However, there are other factors that may possibly af¬ 
fect the redshift of the polarized flux. For example, the broad 
lines are usually redshifted with respect to the narrow emis¬ 
sion lines in many AGNs. The amount of redshift differs for 
various lines, and this may be associated with the stratified 
structures of the BELR consistent with the reverberation 
studies of several Seyfert galaxies (e.g. Peterson 1993). In 
addition to the velocity shift in the source part, the kine¬ 
matics of the scattering agents is important to determine 
the overall profile of the scattered (and polarized) flux. In 
particular, Conway and Blanco (1995) pointed out that the 
HI component in Cyg A is moving with a typical velocity of 
several hundred km s“^. We omit the quantitative analysis 
of this effect in this paper. 

The modelling by Young et al. (1996) indicates that 
the broad lines of IRAS 110548-1131 are polarized by the 
same degree as the continuum, for which case the Raman 
scattering is not required. Therefore, if the Raman scattering 
also operates, then the Raman scattering features may be 
found only in the fine structures in the polarized flux with 
high resolution and good signal to noise ratio requiring large 
amount of integration time. 

Here, one example of the fine structures in the polar¬ 
ized flux can be a relative shift of the narrow line peak and 
that of the polarized flux. The electron-scattered flux re¬ 
duces the asymmetry in the total scattered flux if the broad 
Ha is symmetric. Therefore the relative strength of the elec¬ 
tron (and/or dust)-scattered flux affects the location of the 
polarized flux peak. 

Another possible solution may be provided by observing 
metal lines including Mg H 2800 in the spectropolarimetry 
mode. Because metal lines have negligible scattering cross 
section, they will not be scattered by atomic hydrogen and 
may possess the polarized flux entirely due to free electron. 
Therefore, the Raman scattered Balmer lines can be shown 
by verifying different shapes and widths of the polarized 
fluxes compared with those of metal lines. 

One important process expected from a high column 
neutral hydrogen is the Rayleigh scattering of the photons 
around the Lyman series lines. The Rayleigh scattering of 
Lya photons has been discussed in symbiotic stars by Isliker 



10 Lee and Yun 


et al. (1989). Because the Rayleigh scattering cross section 
is in general larger than the Raman scattering cross section, 
it is expected that the Lyman lines are almost as strongly 
polarized as the Balmer lines. 

In this respect, an interesting observation is the HST 
spectropolarimetry on the quasar PG 1630+337 by Koratkar 
et al. (1995), who found out a strong polarized flux around 
Lya. It is a theoretical possibility that this feature is formed 
by the Rayleigh scattering process, in which case a more 
weakly polarized Ly/I flux is expected. Because the strong 
broad line N V 1240 is present at ~ 6000 km s“^ to the red 
side of Lya, the possibility that the polarized flux is mainly 
contributed from receding N V is not excluded as noted by 
Koratkar et al. (1995) (see also Arav et al. 1995). 

Therefore UV spectropolarimetry on nearby AGN and 
optical spectropolarimetry of high redshift AGN will provide 
strong constraints on the importance of neutral hydrogen 
scatterings. 

5 SUMMARY 

We performed a Monte Carlo computation of the linear po¬ 
larization expected from the Raman scattering process by a 
neutral hydrogen component in active galactic nuclei. The 
main features of the Raman scattering processes in AGN in¬ 
clude a large width due to the enhancement of the Doppler 
shift, and a possible red asymmetric polarized flux around 
Ha ascribed to the asymmetric conversion efficiency with a 
smaller contribution of Ly/I and O vi blended broad emission 
lines. 

It is not still clear whether the neutral component with 
a high column density exists in a typical AGN, and remains 
a very intersting possibility that the Raman scattering pro¬ 
cesses provide a natural interpretation of several spectropo- 
larimetric observations of narrow line AGN including Gyg A 
and IRAS 110548-1131. It is concluded that more detailed 
comparisons of polarized flux of other metal lines and hydro¬ 
gen Balmer lines will put stronger constraints on the possible 
effects of high column neutral components in AGN. 
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